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A sensitive polarization modulation technique uses photoelastic modulation and hetrodyne de- 
tection to simultaneously measure the Faraday rotation and induced ellipticity in light transmitted 
by semiconducting and metallic samples. The frequencies measured are in the mid-infrared and 
correspond to the spectral lines of a CO2 laser. The measured temperature range is continuous and 
extends from 35 to 330K. Measured samples include GaAs and Si substrates, gold and copper films, 
and YBCO and BSCCO high temperature superconductors. 



I. INTRODUCTION 

Conventional dc Hall effect measurements in novel elec- 
tronic materials such as high temperature superconduct- 
ing cuprates (HTSC) have been essential in revealing 
the unusual character of these systems. [Q| Extending 
these measurements to higher frequencies allows one 
to probe more effectively the energy scales of the system 
(e.g., the plasma frequency, the cyclotron frequency and 
the carrier relaxation rates) while minimizing the effects 
of impurity scattering, which can dominate dc measure- 
ments. Further, the Hall angle and Hall conductivity 
satisfy sum rules [? ], which when applied to ac Hall 
data provide insight into the intrinsic electronic struc- 
ture of systems ranging from conventional Fermi liquid 
metals || to more exotic metals such as HTSC and 
other transition metal oxides. || 

Magneto-polarimetry measurements can be used to ex- 
tend Hall effect measurements into the infrared frequency 
range (10 13 Hz). These measurements are sensitive to the 
complex Faraday angle 6*f, which is closely related to the 
complex Hall angle 6*h- || Though 9f and 6*h tend to 
be small for metals in the mid-infrared (MIR, 900-1100 
cm -1 ), there are a number of advantages in performing 
these higher frequency measurements. First, the high fre- 
quency allows one to avoid impurity scattering or grain 
boundary effects which may dominate lower frequency 
Hall measurements. This is especially important in new 
materials which often contain many impurities and de- 
fects. Furthermore, the MIR measurements allow one to 
examine the trends observed at lower frequencies. Since 
tan#H (and $h) obeys a sum rule jT) it is very useful to 



be able to integrate #h to higher frequencies to verify 
whether (and where) the Hall angle sum rule saturates 
or whether there is more relevant physics at even higher 
frequencies. Finally, since the high frequency behavior of 
9u is constrained by the general requirements of response 
functions, a simple, model-independent asymptotic form 
for $h becomes more accurate at higher frequencies. 

In this paper, we present a sensitive mid-infrared po- 
larimetry technique, which uses photo-elastic modulation 
and heterodyne detection, and which can be used to ex- 
plore a variety of materials. || || In what follows we 
shall first introduce the hardware components and de- 
scribe their operation, next we shall present an analysis of 
both the operation and data, and finally we shall present 
recent results on semiconductors, metal films, and high 
temperature superconductors. 



II. EXPERIMENTAL SYSTEM 

A. Overview 

The experimental system of the current work measures 
the very small real and imaginary parts of the complex 
Faraday angle imparted to CO2 laser radiation travel- 
ing perpendicular to and transmitted by the sample film 
which is placed in a perpendicular magnetic field. During 
operation the experimental system performs four major 
tasks essentially simultaneously: 

1) generating and directing a monochromatic, linearly 
polarized light beam normal to the sample 
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2) producing and controlling a magnetic field at the 
sample normal to its surface 

3) setting and maintaining the temperature of the 
sample 

4) analyzing the portion of the light beam transmitted 
by the sample to determine the complex Faraday 
angle. 

Figure |l| schematically illustrates the optical path in 
which the light beam initially issues linearly polarized 
from the CO2 laser and, after enduring various steering 
and attenuation components, proceeds through lens 1, 
which focuses it to a point at the chopper. The chop- 
per impresses onto the intensity of the beam a temporal 
square wave of 112 Hertz which will later facilitate the re- 
moval of laser power variations using ratios. Lens 2 then 
refocuses the emerging beam so as to eventually produce 
a focused diffraction spot on the surface of the sample. 
The intervening element, a film polarizer, "cleans up" the 
beam's polarization, removing any contamination caused 
by the chopper, attenuators, etc. It is this highly linearly 
polarized beam, which after passing through the ZnSe 
magnet housing window encounters the sample as shown 
in Fig. H The sample will reflect, absorb, and trans- 
mit portions of this beam. The portion transmitted will 
have sustained a Faraday rotation proportional to the 
magnetic field and consistent with the physics peculiar 
to the sample. We shall assign the x axis to the initial 
polarization direction and the z axis to the direction of 
propagation. 

The Faraday rotation angle includes both real and 
imaginary terms. The real term corresponds to a simple 
geometric rotation of the polarization vector about the 
direction of propagation. The imaginary term relates di- 
rectly to the ellipticity of the polarization and has been 
called the circular dichroism. The magnitude of the Fara- 
day rotation at CO2 wavelengths is minute — about 10 
radians. As such, one may consider the electric vector 
along the y axis to be a feeble signal added to the com- 
paratively prodigious signal of the original the electric 
vector along the x axis. This small signal cannot be mea- 
sured directly by the common method of cross polarizers 
because the power in the field polarized along the y axis is 
proportional to the square of the Faraday angle or 1CP 8 
times the power in the original beam in the x polariza- 
tion. Clearly, leakage of the input beam through cither 
polarizer would effectively direct power into the y axis 
polarization which would well overwhelm any amount 
occasioned by the sample. The natural choice in such 
instances is heterodyne detection. 

Heterodyne detection produces cross terms which are 
proportional to the electric vector along the y axis rather 
than its square. We realize heterodyne detection by 
phase modulating the electric vector along the y axis 
using a photoelastic modulator (PEM). The PEM modu- 
lator is a model 90 manufactured by Hinds Instruments, 
Inc. with a modulation frequency of 50 kHz. It comprises 



essentially a block of ZnSe and a piezoelectric transducer 
which creates the stress alternations at a frequency of 
50 kHz and ultimately the phase modulation of the light 
beam. The electric vector along the x axis is essentially 
unaffected by the PEM. 

Returning to our discussion of the optical path the 
beam emerges from the PEM and then strikes a polar- 
izer on a 45° angle, which, simply stated, allows a portion 
of the electric vector from each of the y and x directions 
to contribute to that which emerges now at a 45° an- 
gle. These two contributions after being mixed by the 
detector produce the 50 kHz PEM frequency along with 
sidebands. The amplitude of each of the sidebands is pro- 
portional to the electric vector along the y axis and also 
proportional to the corresponding Bessel function whose 
argument is the depth of modulation. In fact, as will 
be derived below, the even sidebands or harmonics are 
proportional to the real part of the Faraday angle and 
the odd sidebands or harmonics are proportional to the 
imaginary part of the Faraday angle. 

In the following sections we consider certain elements 
or subassemblies of the experimental system in detail be- 
ginning with the CO2 laser. 



B. C0 2 Laser 

The CO2 laser, fashioned from a model PL5 manu- 
factured by Edinburgh, provides a number of spectral 
lines ranging from 9.174 to 10.860 microns (920 to 1090 
cm -1 ), any one of which is selected by a grating internal 
to the cavity. A separate CO2 laser spectrometer ver- 
ifies the wavelength of each line. The direction of the 
beam exiting the laser, however, differs for different lines 
sometimes by nearly as much as one half of a spatial 
mode. Without realignment the result would be a spa- 
tial shift of the diffraction spot at the sample by as much 
as one half of a spot diameter. The small size of a sample 
can aggravate the sensitivity of the system to such slight 
misalignment particularly when such changes cause the 
beam to wander off or even near the edge of the sam- 
ple. Lens 3 essentially focuses an image of the sample 
onto the iris shown in Fig. [l] and may be used as and aid 
in regaining alignment as is necessary to compile valid 
relative frequency data. The design and construction of 
the hardware well ensures that the relative position of 
the sample, lens 3, and the iris are constant. Therefore, 
adjusting the laser beam steering to maximize the power 
through the iris after each laser line change likewise en- 
sures that the beam is passing through the same area of 
the sample. While actually taking data the iris is dilated 
to accommodate a small amount of beam wander. 

Another consideration of the laser is the high output 
power — 2 to 30 Watts. Three attenuators reduce the 
power ~30dB to avoid heating the sample while still pro- 
ducing a robust signal within the linear range of the MCT 
detector. These attenuators are ZnSe windows having an 
antireflection coating on one side. They replace mirrors 



3 



in the beam steering assemblies with the uncoated side 
serving as the reflective surface. Each of these attenua- 
tors introduces approximately 10 dB of loss. 

The output power of the laser, notwithstanding its har- 
diness, fluctuates temporally. Therefore, producing us- 
able data requires the formation of simultaneous ratios 
between the 2nd and 3rd harmonics, and a sample of 
the laser power level. The current experiment employs 
the usual method of "source compensation" |Q which 
involves the chopper along with an additional lockin am- 
plifier. Note also that the heterodyne detection system 
described above permits use of the same detector for the 
power level sample as well as the 2nd and 3rd harmonics. 
As to be derived below, the power level sample, aside 
from the chopper frequency, corresponds to the original 
optical carrier frequency centered within the sidebands 
mentioned above. This feature eliminates the wavelength 
dependence of the detector and thus facilitates accurate 
relative measurements of the wavelength dependence of 
the complex Faraday angle. A final consideration regard- 
ing the laser involves its placement with respect to the 
magnet. The particular position of the laser is perpen- 
dicular to the magnetic field and at such a distance as 
to reduce detuning of the laser cavity caused by magne- 
tostriction to a tolerable amount. 



C. Magnet System and Sample Mounting 

The magnet system began as an 8 Tesla, split coil, 
Helmholtz Spectromag manufactured by Oxford Instru- 
ments and, except for modifications to the external win- 
dows, the sample handling hardware, and the internal 
bore tube shielding, it has served the experiment with lit- 
tle modification. To understand the modifications to the 
external windows, consider that the wavelength range of 
interest requires ZnSe windows. However, the band gap 
of ZnSe is only 2.7eV, and the magnetic field of 2 Tesla at 
the original location of the windows caused a significant 
Faraday rotation, which overwhelmed that of the sample. 
The extension tubes locate the windows where the mag- 
netic field is less than 0.1 Tesla. The remaining Faraday 
rotation background is thus reduced to the same order as 
that of the sample and having been carefully measured 
for different wavelengths can be easily removed from the 
data at a later time. For magnets requiring cold windows 
located within the high field region the background may 
be reduced by using a material with a higher band gap 
such as BaF. The sample handling hardware, ill-suited 
to the task at hand, enjoyed substantial redesign and re- 
manufacture. The effort included adding a steady pin 
and vise, implementing a more responsive temperature 
control system, and developing a near stress free sam- 
ple mount. The steady pin, depicted in Fig. || protrudes 
from the blade of the original sample stick to which it 
is brazed. The vise, located within the magnet hous- 
ing, engages the pin using a teflon collet. This combined 
apparatus restrains the sample stick against the forces in- 



duced by the magnetic field and thus prevents the large, 
extraneous interference signals resulting from changing 
multi-paths and etalons. 

Figure || also depicts parts of the new temperature 
control system which comprises a heater and a cooling 
link. The heater is a 400 Ohm 1/4 Watt metal film resis- 
tor potted into the copper sample carrier using 2850FT 
Epoxy with Catalyst 9 both of which are manufactured 
by Emerson & Cuming. The cooling link is simply a 4 
cm length of 22 gauge copper wire connecting the copper 
sample carrier to the blade of the original sample stick 
from which it is otherwise thermally insulated by spacers 
and nylon screws. The operation is delightfully simple: 
Liquid helium delivered to the sample stick by the orig- 
inal provision, cools the blade to about 10 K, and the 
cooling link cools the copper sample carrier. Current de- 
livered to the resistor in an easily controlled fashion can 
provide up to 2 Watts of heat to the copper sample car- 
rier. Because of the low thermal mass of both the copper 
sample carrier and the link, a compromise temperature 
emerges within about one minute. Sweeping the temper- 
ature entails nothing more than adjusting the current. 
A important advantage of this system is that only the 
rather small copper sample carrier changes temperature. 
The remaining hardware, remaining essentially constant 
in temperature, finds little reason to move or warp so as 
to adversely affect the measurement. 

The near stress-free sample mount consists of a phos- 
phor bronze wire retaining spring and a thermally con- 
ductive flexible silver-filled RTV known as Eccobond 
59C manufactured by Emerson & Cuming, Inc. The 
particular sample shown attached to the sample stick 
is a small irregularly shaped film of optimally-doped 
Bi2Sr2CaiCu208+5 approximately 200 nanometers thick 
which had been peeled from a single bulk crystal. This 
small film was placed against a polished surface of a bar- 
ium fluoride crystal which serves as a substrate and main- 
tains the temperature of the sample. Van der Waal's 
force holds the film in place. A 2° wedge of the BaF 
substrate eliminates etalon effects. Only one of the two 
corners of the substrate is cemented to the copper sam- 
ple carrier. The phosphor bronze retaining spring lightly 
holds the other corner while allowing some motion to re- 
lieve the stress caused by dissimilar thermal expansion 
coefficients. Without this provision the stress induced 
in the substrate had caused overwhelming and unpre- 
dictable complex Faraday rotations. It is this stress free 
mounting which actually facilitated the fast temperature 
scans required to eliminate 1/f noise apparent in previous 
work. The copper sample carrier is provided with a small 
indexing hole also shown in Fig. [j[ Prior to insertion of 
the sample stick into the magnet the sample's position 
is accurately measured with respect to this hole. After 
insertion, the indexing hole is located using the trans- 
mitted intensity of the laser beam. The sample is then 
positioned within the beam by raising or lowering the 
magnet and sliding the sample stick in or out. Finally 
the new reentrant bore tubes are fit with graphite plates 
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to absorb the stray radiation scattered from the incident 
laser beam by the various reflective surfaces, e.g., the 
sample, the substrate, the magnet windows, etc. 



D. Optical Table Components 

The remaining items of Fig. [j] which along with the 
magnet are located on the aluminum optical table will 
receive consideration within this section. Among these 
items are the chopper and lens 1. As explained above, 
the chopper impresses upon the beam a square wave am- 
plitude variation of 112 Hertz. Lens 1 participates in 
this task by focusing the beam to a point at the chopper 
blades. Thus, the passage of a blade cuts the beam off 
and on in an most abrupt manner. This prevents spa- 
tial variations within the beam from causing phase and 
amplitude errors in the reported laser power level. 

Also of particular importance a pivot platform upon 
which sits the optical detection system. Two aluminum 
box beams connect this platform to a vertical pivot lo- 
cated directly below the sample. When the two clamps 
securing the platform to the optical table are loosened 
the platform may be rotated about the vertical axis of 
the pivot. This motion is necessary to align the detection 
system to both the indexing hole and wedged substrates 
which bend the laser beam about the same axis by an 
amount depending on the substrate wedge angle and its 
index of refraction. 



E. Optical Detection System 

The first component encountered by the laser beam af- 
ter passing through the magnet window is lens 3 which 
focuses an image of the sample onto the iris as previ- 
ously discussed. This lens joins with the PEM and po- 
larization analyzer in an assembly, which can be rotated 
as a single unit about the input optical axis, which is 
the z axis. With a moment's thought one will recognize 
that such a rotation is equivalent to a real, but opposite, 
Faraday rotation, at least for sufficiently small angles so 
that the polarization sensitivity of the MCT detector is 
not apparent. Fixed calibration rotation stops limit this 
rotation to a known amount thus serving as a Faraday 
rotation calibrator. Simply put, one rotates the assembly 
the known amount and uses this to scale the empirical 
values for each wavelength. 

Another important consideration involves reflection at 
the surfaces of the ZnSe interaction block of the PEM. 
A reflected beam which makes additional passes through 
the ZnSe interaction block will receive additional modu- 
lation. Since, the cross term is a function of the depth 
of modulation, such triply modulated stray beams can 
cause significant errors. An AR coating and a tilting of 
the PEM by 25° reduce and displace reflected beams and, 
thereby, sufficiently reduce their effect. 



A variable selection of polyimide films further attenu- 
ate the laser beam so that its power is within the linear 
range of the MCT detector which is a model J15D14 mer- 
cury cadmium telluride detector manufactured by EG&G 
Judson. 



F. Electronic Instrumentation 

Along with the electronics which attend the mag- 
net, PEM, chopper, etc., the system uses three model 
7260 harmonic lockin amplifier/detectors manufactured 
by EG&G. These lockins differ from earlier styles in that 
they detect and report the RMS voltage of a selectable 
harmonic of the input signal. This feature is essential 
because, as mentioned, the even harmonics of the phase 
modulation of the PEM are proportional to the real part 
of the Faraday angle and the odd harmonics are propor- 
tional to the imaginary part. The first lockin determines 
the RMS voltage at 112 Hertz and is usually considered 
to be the DC reference. The second lockin determines 
the RMS voltage at 2u>, the 2nd harmonic of the PEM 
frequency, and the third determines the RMS voltage at 
3u>, the 3rd harmonic. A Labview program running on 
a local computer records data points consisting of the 
following: 

t> n ii V 2u (x) V 2ul {y) 
time or B field, , , 

Vll2 Vll2 



Vi 



-, V u 



112 



If the beam moves within the PEM aperture the phase 
of the harmonics will change. This phase is initially set 
such that the y channel is very nearly zero with most of 
the signal in the x channel. To avoid any errors from the 
changing phase we use the positive definite magnitude 
formed from the combination of the x and y channels 
and apply the sign chosen determined by the calibration 
procedure described later. 



III. ANALYSIS 

In the analysis of both the system and data we employ 
two formalisms. One concerns the representation and 
transformation of the polarization state of light and the 
other concerns the transmission and reflection response 
of multilayer stacks with complex material properties, cr, 
e, or fi. Both are presented as appendices. 



A. Faraday Angle 6f versus Lockin Outputs 

For light initially propagating in the z direction and 
polarized along the x axis, we define the Faraday angle 
as 



arctan ■ 



txx 



v yx 

txx 



(1) 
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where t xx is the field transmission along the x polariza- 
tion, t yx is the field transmission along the y polarization, 
and the arctan function was dropped because of the very 
small angles involved in the current work. As such it rep- 
resents the complex amplitude (amplitude and phase) in 
the y polarization having been derived from the incident 
radiation in the x channel. One will recognize Re (Op) as 
the geometric rotation of the incoming polarization and 
Im (Of) as the circular dichroism. 

Consider next that the experimental system examines 
the light transmitted by the sample, which results from 
incident light initially polarized along the x axis (recall 
the z axis is the direction of propagation). A simple 
matrix equation represents this activity in the linear po- 
larization basis as 







•^out 
2/out 



(2) 



Assuming the a and b axes of the sample to be indistin- 
guishable, i.e., near square symmetry, and noting that 
the B field is uniform and parallel to the z axis, we know 
that the sample transmission of Eq. (|J) will be diagonal 
in the circular basis. Using the polarization formalism in 
the appendix the sample can be represented as 



t p 

o t n 



t, 



% to. 









H - ^ tyx 



(3) 



where t p is the transmission coefficient corresponding to 

positive helicity (positive rotation seen at fixed point 
for a wave traveling in the positive z direction, and t n 
is the transmission coefficient corresponding to negative 
helicity (negative <fi rotation seen at fixed point for a wave 
also traveling in the positive z direction. Transforming 
back to the linear basis the sample matrix becomes 



/ txx 


t X y \ 




( txX tyx \ 




tyy ) 




{tyx txx ) 



which is considerably simpler. It is also useful to note 
that by the definition in Eq. ([!]) the Faraday angle relates 
very simply to the ratio of the transmissions of the right 
and left circular polarization: 
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Next, for the purpose of developing a relation between the 
sample transmission and the lockin outputs consider the 
relevant experimental elements represented schematically 
as: 



laser — * sample — * PEM 
— > square law detector - 



polarizer at 45° 
lockins. 



Again using the polarization formalism, the components 
in the the first line, which operate on the laser beam |x) 
may be represented as: 



R © 



1 




x R 



— ) x PEM x sample x \x) 



(p|out) 
(nlout) 



~ l - iflOWe'i 

J y J \ e- i2 i 

1 -i \ ( e if3cosuJt \ / 1 1 



1 i 



1 



—lap Q 

e l8F 



i -n / 1 

1 i \ 



where ui = 2ttJ is the radial frequency of the PEM modu- 
lation and (3 is the PEM modulation depth or amplitude. 
The above signal is incident upon a square-law detector 
whose output is a voltage proportional to the square of 
the modulus of the amplitude of the fields: 

voltage oc (out | out) oc optical power. 

For the small angles considered the multiplication of 
these matrices along with the Bessel function expansion 
of e i/3cos(27r/) produces 



Re( 



Im ( 



Re 



Im 



v yx 

txx 
tyx 

txx 



1 



2 a.' 



4J 2 (/3) V dc 

1 Vzu> 
iJ 3 (p) Vd c 



(4) 
(5) 



where Ja(/3) and Js(f3) refer to 2nd and 3rd order Bessel 
functions respectively, Vdc refers to the optical power at 
dc (actually at 112 Hz due to the chopper). 



B. Sample Properties versus Complex Faraday 
Angle 

Since the sample transmission matrix, Eq. (||) is di- 
agonal in circular polarization, each circular polarization 
channel acts independently and may be treated as such. 
Indeed, by the same symmetry considerations, the sam- 
ple matrix representation of the conductivity tensor is 
also diagonal in the circular basis. In what follows we will 
develop an expression for the transmission of a sample- 
substrate combination and then separately apply this ex- 
pression to each polarization channel. This will give us 
an expression for the indices of refraction of the film for 
positive and negative helicity rif p and nf n respectively, 
which are then easily related to a xx and a xy . 

Figure ^| appearing earlier depicts the transmission 
through the sample, which is usually a film-substrate 
combination. The input beam strikes the sample film, 
a portion reflects as shown, and a portion propagates 
into the sample film and after enduring some absorp- 
tion arrives in the forward direction at the output side of 
the film. This portion strikes the interface between the 
film and the substrate mounting surface. Some reflects 
back and forth within the sample film and some proceeds 
as a beam into the wedged substrate soon reaching the 
substrate-air interface. Most of this beam propagates 
into the air eventually arriving at the detection system. 
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Because of the wedge angle of the substrate-air inter- 
face, the small reflected portion leaves the optical path 
and is absorbed by strategically placed graphite slabs lo- 
cated in the reentrant bore tubes. For this reason the 
substrate-air interface does not participate in the Fara- 
day rotation for any of the materials of current interest 
and need not be considered. All of the relevant activity 
above can therefore be represented schematically as 

air . film . substrate 



The formalism for multilayer transmission presented in 
the appendix represents this sequence generally as 



(sample) 



E n 



( -E'sub, 

I 



(6) 



where (sample) = S fi i m , sub = UfiimS a ir,fiim with S denot- 
ing the interface matrix and U denoting the propagation 
matrix. That is, the sample acts on the electric field 
vectors and _E a j r j4 _ of the incoming and reflected 

waves respectively on the input side producing the elec- 
tric field vectors £ su b,— » and £ su b,<- of the outgoing and 
reflected waves on the output side, where, of course, we 
have no incoming wave. We desire only the transmission 
which from Eq. (||) is the (1,1) element of the inverse of 
the sample matrix. Since the transmission is the same 
for either propagation direction; a moments thought will 
lead one to propagate backwards to avoid taking the in- 
verse of the sample matrix. The end result is the trans- 
mission is given by the inverse of the (2, 2) element of the 
sample matrix, which using the relations in the appendix 
can be related (after some algebraic manipulation) to the 
material properties: 

t(n) = An s n f e ikd 

i2kd 



x [n s (e rzkd (n f -l)+n f + l) 
-n f (e^Mfo-iJ-n/-!)]- 1 



(7) 



where d is the thickness of the film, k — 2im / /A is the 
wavenumber within the film, and n/ and n s are the in- 
dices of refraction of the film and substrate. There is 
one such equation for each circular polarization channel 
n and p. 

Analysis of the data ultimately requires an equation for 
the index of refraction (or some other material property) 
in terms of the Faraday angle. Equation (J3J) with Eq. (Q) 
contain the pertinent information but the combination 
does not lend itself to inversion. Fortunately, because 
the relative difference between the indices of refraction of 
the film for n and p polarizations is usually very small, 
we can generate a readily invertible form by expanding 
their combination about either nf >n or Uf tP . Choosing 
n f,p ~ n f w ith S — nf, n — n>f tP , expanding Eq, (||), and 
keeping only the linear term results in 



-id 1 dt(n) 
2 t(rif) dn 



(8) 



Combining Eqs. (^) and (@) produces 



26» F 



r?j (n s + 1) —inf tan (kd) ^ 



n s + n' 



kd 



(n s + n^J — tan (kd) (n s + ikd (n s + 1) nf — rij\ 



(9) 



For an expression for a xx and <j xy note that by the same 
symmetry assumed earlier the film dielectric permitiv- 
ity tensor and the conductivity tensor are also diagonal 
in circular polarization, which we can convert to linear 
polarization as 



(C\a\C) 



(C\L 




L){L\C) 



Using Maxwell's equations the complex conductivity in 
CCS units is therefore 



ff «=-ft?( B /,« 



<Jxx = - 3- ( r 



■f,n 



Ait 



x ^ 
= O—Uf. 

47T 



(10) 



(11) 



By means of Eqs. © and ©, Eq. (§ becomes 



>xy 



2tt 



if (n s + 1) — i tan(fcef) (^n s + n 2 ^ 



kd (n s + nj^j — tan(fcd) (n s + ikd (1 + n s ) nf 



(12) 



This expression of course requires accurate values for the 
real and imaginary parts of Q-p. The experimental appa- 
ratus, however, introduces several errors into these val- 
ues, which we have ignored up to this point but, which 
nevertheless must be removed or compensated as de- 
scribed below. 



C. Calibration 

There are three rather serious sources of error ignored 
in the foregoing. First, the MCT detector and subsequent 
electronics introduce their own response functions which 
must be removed from the data. Second, the PEM mod- 
ulation reported by the PEM electronics is frequently in 
error for reasons to be discussed below. Third, the Fara- 
day angle 9f determined by the above equations is the 
total Faraday angle, which includes a background, which 
results mostly from the ZnSe magnet windows and the 
substrate. 
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Regarding the response functions, Eqs. (Q) and (||) 
above require the ratio of the optical power in the 2nd 
and 3rd harmonics to the optical power at 112 Hertz 
(referred to as DC) at the MCT detector. However, 
the MCT and other electronics possess frequency trans- 
fer functions, which attenuate these signals by different 
amounts, and once determined must be divided out of the 
collected data. These functions were determined using a 
high-speed communications type laser diode to generate 
a known optical beam at the same frequency and ampli- 
tude as that within the experimental system. When this 
beam was directed at the MCT the result corresponded 
rather well to the response of a simple RC circuit. 

Next consider that the first lockin measures the RMS 
voltage of only the 1st harmonic of the 112 Hz square 
wave chopped signal, whereas Eqs. (Q) and (||) require 
the dc power. To determine this correction factor we 
examine the relation between a square wave and its first 
harmonic. The normalization requirement generates the 
inverse of the first harmonic: 



r\t)f(t) dt = i 



(13) 



Therefore, representing the harmonic as sin(urf), the nor- 
malization constant N is determined by: 



f~\t) sin(wi) dt = N r {t) sin(urt) dt = 1 
/o Jo 

(14) 

Therefore N — - and / 1 (t) — -sm(uit). We now per- 
form the integral transform to determine the first har- 
monic amplitude for a square wave of peak amplitude 
1. 



/ f (*) (square wave) dt = — 
Jo * 



(15) 



So, for a square wave of peak amplitude 1, the output 
of the 7265 would be the RMS value for a sinewave of 
peak amplitude 4/tt. Because data2cj and data3w (also 
recorded in RMS) were divided by the 112 Hz RMS 
signal, removing this effect therefore entails multiplying 
them by 4/-7T. 

Next, consider the PEM. The PEM electronics as sup- 
plied by the manufacturer at best only infer the ampli- 
tude of the phase modulation impressed upon the trans- 
mitted beam. This inference does not consider the path 
of the beam through the ZnSe interaction crystal, the 
ZnSe crystal temperature, nor the possible changes in 
coupling between the driving quartz crystals and the 
ZnSe interaction crystal. Equations (Q) and (^) , how- 
ever, require an accurate knowledge of this modulation 
in the value of the variable fj. The current experiment 
explores, among other things, the frequency dependence 
of 9p. Because the frequency dependence is often weak, 
its determination requires data with very high relative 
accuracy. A simple method arises for obtaining Re (6-p) 
which makes use of the provision for physically rotat- 
ing the PEM assembly back and forth a predetermined 



amount as discussed above. Prior to collecting each data 
set for each laser line one rotates the PEM assembly back 
and forth by the predetermined amount and records the 
readings from the lockin amplifiers. This defines the cal- 
ibration factor for Re (Op): 



C 



real 



predetermined physical rotation 
reading2 — readingi 



where 



readings = 



readmg(2w) 



(16) 



(17) 



reading(112)i 

This calibration procedure also supplies the sign of 
Re (#f), which the electronics by itself leaves somewhat 
ambiguous. Subsequent to this calibration the real part 
of the Faraday angle emerges from the data as 



Re( 



(sign) x data2w x C r 



cal • 



(18) 



Because the frequency dependence of Im is also 
weak, its determination also requires data with very high 
relative accuracy. Equation (|4|) along with the rotation 
calibration above also afford a relative calibration for 9^. 
This calibration is essential in order to reliably deter- 
mine the correspondingly small wavelength dependence 
of Im (9-p) The procedure entails simply comparing the 
actual, predetermined, physical rotation with the calcu- 
lated value from Eq. (Q) and then adjusting the value 
of (3, the retardance until they agree. When calculat- 
ing Im(0 F ) with Eq. (§) we use this adjusted value /3 a dj 
in place of that determined from the manufacturers cal- 
ibration. This procedure should be performed initially, 
and then after each laser line change. Here again the 
electronics leave the sign ambiguous, so it must be de- 
termined using a sample of known ellipticity such as a 
quarter waveplate. Combining the foregoing calibration 
corrections: 



Im( 



= (sign) 



(/3adj) 



data3w(A). (19) 



Finally consider the background contribution to 9f 
which manifests itself as a number of terms: 



'F, background 

= cq + c x B + c 2 B 



(20) 



where each coefficient Cj may have a wavelength depen- 
dence. In the current work the data sets contain 9f as a 
function of magnetic field which is, for example, scanned 
from +8 to -8 Tesla or vice versa. From this data we com- 
pute the slope 89f / dB which obviously does not contain 
c . Further, because the magnetic field values are both 
positive and negative, computation of the average slope 
eliminates all terms even in B. Removing the remain- 
ing odd terms requires direct measurement of the back- 
ground using a sample of the substrate material followed 
by simple subtraction of the result from the data sets. 
The results of these measurements for substrates of Baf 
and high purity Si were purely real Faraday angles, which 
were consistent with the bandgap of these substrate ma- 
terials and the ZnSe magnet housing windows as well as 
the absence of free carriers. 
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IV. EXPERIMENTAL MEASUREMENTS 
A. Semiconductors 

Faraday measurements were first performed on semi- 
conductor samples for three reasons. First, the signals 
are large and can be readily increased by using thicker 
samples. Second, the optical properties of semiconduc- 
tors are well known, which allows the accuracy of the 
Faraday measurements to be verified by by comparing 
them to the results in other experiments. Finally, since 
the thin film samples are grown on semiconductors, it is 
important to measure their contribution to the Faraday 
signals accurately in order to remove this background 
from the desired thin-film signal. 

The Faraday measurements were performed on three 
semiconductors. In semiconductors, two sources con- 
tribute to the Faraday signals plf . The first source is 
from free carriers, which have a contribution to the Fara- 
day signals that is proportional to w~ 2 , where ui is the 
radiation frequency. The second source is from interband 
transitions, where the magnetic field causes anisotropy in 
the refractive indices for left and right circularly polar- 
ized light. This difference in index leads to a phase shift 
between left and right circularly polarized light, which 
in turn results in a rotation of the incident polarization 
(Faraday rotation). For below band gap radiation, the 
interband contribution is proportional to uj 2 . Since the 
MIR frequencies are an order of magnitude smaller than 
the energy band gap in typical semiconductors, no inter- 
band absorption occurs and therefore only the interband 
term contributes to Faraday rotation. Furthermore, since 
the semiconductors in this experiment are of high purity 
and the MIR frequencies are relatively high, the free car- 
rier contribution is minimized. Faraday rotation caused 
by such samples is quantified in terms of its Verdet con- 
stant V, which is defined as the angle of rotation per 
unit magnetic field per unit thickness of sample. The 
frequency dependence of the Verdet constant and its rela- 
tionship to the Faraday angle at 8 T for a semiconductor 
is given by the following equation []I2|: 



V 



ULO 



v _ Re[0 F (8 T)] 180 
^2 ~~ 8TxD X ~ 



(21) 



where u is the coefficient for the interband contribution 
and v is the coefficient for the free carrier contribution 
to the Faraday rotation, and D is the sample thickness. 
Rcf. [l2] et al. have measured the Verdet constant as a 
function of frequency (2000-3300 cm -1 for several Si sam- 
ples with different free carrier concentrations. One can 
extend these results to 949 cm -1 by using Eq. ^l] to ob- 
tain u and v from the data in Ref. [lj. The Verdet coeffi- 
cient at 949 cm -1 is estimated to be 25.0 degree/T-m. 
The value obtained at that frequency by MIR Fara- 
day measurements on a high purity (> 1 kfJcm) sil- 
icon sample is 25.6 degree/T-m. This corresponds to 
Rc[0 F ] = 2 x 10" 3 at 8 T and 949 cm^ 1 for a 0.50 mm 
thick sample 



The second semiconductor material is semi-insulating 
GaAs. In this case, the frequency dependence of the 
Faraday rotation signal is explored. Fig. [| plots the 
Re[#p] as a function of frequency squared at 8 T for a 
GaAs sample. The solid line represents a w 2 -fit. The 
semi-insulating sample is wedged with an average thick- 
ness of 0.42 mm. The ui 2 frequency behavior of Re[#F] is 
consistent with Eq. |l|. At 8 T and 949 cm" 1 , the Re[0 F ] 
is approximately 3 x 10~ 3 . The Verdet coefficient is ap- 
proximately 44 degree/T-m. The value obtained at 10.59 
fim for Cr-doped|l3j GaAs by Ref. [l4] is 41 degree/T-m. 
This is approximately 50 % larger than that for Si, which 
is consistent with the fact that GaAs has direct band gap 
that is approximately a factor of two smaller than that 
of Si. This brings the MIR radiation in GaAs closer to 
the refractive index anisotropy at the band edge. In fact, 
the Re[6>F] for GaAs decreased as temperature is lowered, 
which is consistent with the increase in the band gap as 
temperature decreases. As expected, no circular dichro- 
ism signal is observed (Iiti^f] ~ 0). 

The final semiconductor investigated is LaSrGaO.4. 
This material has the largest band gap, and hence the 
smallest interband contribution to Re[0 F ]- At 8 T and 
949 cm -1 , the Re[6*p] is approximately 6.5 x 10~ 5 for a 
0.31 mm thick sample. This translates in a Verdet coef- 
ficient at 949 cm -1 of 1.5 degree/T-m. Despite the onset 
of strong absorption by phonons for radiation below 1000 
cm^ 1 , LaSrGa04 proved to be an excellent substrate due 
to its weak contribution to the Re[0F] background. Note 
that as frequency is decreased, the strength of the semi- 
conductor Faraday signals decreases (Eq. ^l]) while the 
signal from the free carriers in the thin-film metallic sam- 
ples increases (Eq. ||l]). As a result, the lower frequency 
measurements are simpler and more accurate. 



B. High Temperature Superconductor 

Bi 2 Sr2CaiCu208+i 

The measurements also include the response of 
Bi 2 Sr 2 CaiCu208+5 to input radiation from 920 to 1090 
cm -1 and over a temperature range from 30 to 330 K in 
an external B field ranging from -8 to +8 Tesla. Figure ^ 
depicts result 9-p versus temperature. This response is 
not too unlike the Hall angle #h = &xyl&xx with which 
it may be compared (see for example Fig. |?]). To form 
#h one would usually divide a xy by the value of a xx cor- 
responding to the same temperature and frequency. But 
here, we wish to examine that part of the Hall angle re- 
lated to the "free carriers." Interband transitions which 
are on the order of 1 volt contribute essentially nothing 
to a xy . However, we must remove their somewhat more 
substantial contribution to a xx before taking the ratio. 
In CGS units 



Ctotal — Abound + 4" 



. 47T(7f r . 



.47TC7- 



total 



CO 



(J 



(22) 
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<7fr 



"total 



4-7T 



(23) 



Figure ^ displays the values of #h versus temperature. 
The real part of #h increases with temperature but ap- 
pears to saturate around 300 K. Note that the real part 
of the Hall angle is positive for all temperatures in the 
normal state and does not show any striking discontinu- 
ity at T c at 90 K. This data appears in Ref ^ where it 
is compared with far-infrared data from Ref [13. The 
results show a significant disconnect from the behavior 
of the Hall angle in the existing data for YBCO in the 
far-infrared, which indicate a negative value for the real 
part of the Hall angle above 250 cm -1 , whereas that of 
the current work at 1000 cm -1 is positive. The current 
work when analyzed using an extended Drude formalism 
results in a Hall mass comparable to the ARPES Fermi 
mass and a scattering rate comparable to the DC longi- 
tudinal, DC Hall, and far-infrared Hall scattering rates, 
which, however, are only 1/4 of the ARPES values. 



a relaxation rate behavior based on quasi-elastic scatter- 
ing from superconducting fluctuations in the normal state 
of high T c materials. Fluctuation effects have also been 
observed in normal state of underdoped Bi2Si'2CaCu20s 
in measurements of the longitudinal conductivity by THz 
spectroscopy. |2^] The frequency and temperature depen- 
dence of the quasiparticle Hall response serves as a crit- 
ical test for any theoretical model that attempts to ex- 
plain the nature of the normal state in HTSC. 
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C. Metal Films 

Figure ^ shows the 2uj m and 3w m normalized signals 
as a function of magnetic field B at room temperature 
for a Cu thin-film sample. The MIR radiation frequency 
is 949 cm" 1 . The background contribution due to the 
substrate and windows has not been removed. Note the 
negative offset in the 2u> m data. The signals are linear in 
B, as expected. 

It is interesting to compare the MIR Hall effect results 
on HTSC|, | [U| with those on Au and Cu films |, g. 
Surprisingly, the Hall angle results for both materials are 
qualitatively similar, and are well parameterized by the 
Drude form for 6*h- For both systems, the scattering rate 
7h associated with the Hall angle is linear in temperature 
and independent of frequency. In Au and Cu the longi- 
tudinal scattering rate "f x x (obtained from infrared con- 
ductivity measurements) is also temperature dependent 
and frequency independent as expected from electron- 
phonon scattering when the measurement frequency is 
higher than the Debye frequency. For optimally doped 
YBa2Cu307 in the MIR, however, while 7h has a strong 
temperature dependence and no frequency dependence, 
7a; x is temperature independent but frequency depen- 
dent (as found in infrared conductivi ty jl7[ and angular- 
resolved photoemission spectroscopy |18| ) The behavior 
of j x x precludes phonons or magnons as the dominant 
scatterers and the lack of frequency dependence of 7h is 
in contrast to the predicted and observed behavior of a 
Fermi liquid or inelastic scattering in general. Therefore 
the frequency and temperature dependence of 7h that 
are reported in the MIR^] are highly unusual and in- 
dicate a non-Fermi liquid behavior of the normal state 
of YBa2Cu307. Similar results have been obtained re- 
cently in other HTSC materials such as optimally doped 
Bi 2 Sr 2 CaCu20 8 § and underdoped YBa 2 Cu 3 07.[[l6) We 
note that Ioffe and Millis [jl9| have recently proposed such 



APPENDIX A: POLARIZATION ANALYSIS 
USING DIRAC NOTATION 

In polarimetric measurements a train of optical com- 
ponents, along with a sample, modify or operate on an 
input beam, which is in some state of known intensity and 
polarization. Optical components typically include some- 
what imperfect devices such as polarizers, waveplates, 
photoelastic modulators, and mechanical rotators. The 
modified output beam generally strikes a detector pro- 
ducing an output electrical signal. Here we develop a 
general formalism useful in relating the electrical signal 
to the properties of the sample. The formalism uses Dirac 
notation which well differentiates between a vector or 
state \tp) and its representation in some basis: (x\ip). In 
polarimetry the usual bases are linear polarization and 
circular polarization. The formalism also replaces the 
Jones matrices with generalized operators and introduces 
the basis transformations and geometrical rotation oper- 
ations. Finally the result is related to the output of the 
typical optical, "square law" detector. 

In what follows the direction of propagation is in the 
positive z direction and all angles are measured in the 
positive radial direction off the x axis and about the z 
axis. Beginning with the input beam, the formalism sim- 
ply considers it to be a non-normalized ket 

|in). 

Expressed in the linear polarization basis it is 
|in)=^|L)(L|m) = |x)( : r|m) + |y)(y|m) 

L 

or in the circular basis 

|in) = ^|C)<C|m) = b)(p|in) + |n)(n|m) 
c 
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where \p) indicates an electric vector rotating in the pos- 
itive radial direction about the z axis as seen at a fixed 
point and |n) indicates an electric vector rotating in the 
negative radial direction about the z axis also as seen at 
a fixed point. \p) and \n) are also referred to as having 
positive and negative hclicity. In matrix notation 



(L|in) 



E,, 



and 



<C|in> 



The various components in the optical train appear as 
operators acting on the input ket to produce the output 
ket: 

|out) =TVW|in). 

In the linear basis this is 

|out) = \L){L\TVW\L'\m) 

= \L) (L\T\L") (L"\V\L"') (L"'\W\L'} (L'\m) 
= \L)(L\T\L)(L\V\L)(L\W\L)(L\m) 

where summation is assumed over repeated state designa- 
tions L; and the primes have been dropped with the un- 
derstanding the order of the matrices will not be changed. 

Some devices or operations such as geometric rotation 
and Faraday rotation are most easily represented in the 
circular basis. The transformation matrices, (C\L) and 
(L\C) , are the means of conversion between these two 
bases. The transformation for states proceeds as 

|in)= 5>>< L R 

L 

= -£\Q(C\L)(L\m) 

LC 

= \C)(C\L)(L\m), 



and for operators as 

T= ^2\L)(L\T\L')(L'\ 



L,L' 



= S 2 \ C )( C \L)(L\T\L')(L'\C)(C'\ 

L,L' C,C 

= \C)(C\L)(L\T\L)(L\C)(C\ 



In matrix representation 



(C\L) 



V2 U 



and (L\C) 



J_/l 1 

V2 V* -* 



Below are the representations of some simple optical com- 
ponents in the linear bases: a perfect polarizer aligned 
along the x axis 

1 


a 1/4 waveplate aligned with slow axis along x axis 



e 2 




a photoelastic modulator aligned along the x axis 

e i/3cos(27r/t) Q 







1 



Another concern for the devices or components is their 
orientation. Assume, for example, that the operator T 
represents the effect of some component such a polarizer 
which is aligned with its critical direction at an angle 6 
from the x axis. A simple rotation operation develops an 
expression for the rotated polarizer: 

T{6) = RTR 1 

where R is an active rotation and in circular representa- 
tion is 



(C\K(6)\C) 



e~ ie 
e lf> 



Finally, the output of a square-law type optical detec- 
tor, aside from any responsivity factor is for our example 

d(t) cx (out | out) 

where \out) is the Hermitian conjugate of {out\. 



APPENDIX B: MULTILAYER RESPONSE USING 
RELATIVE IMPEDANCE MATRICES 

Assume that some layered materials are arranged nor- 
mal to the z axis and we desire a formalism for the trans- 
mitted and reflected intensities for plane waves incident 
at arbitrary angles. To begin, the plane wave solution to 
Maxwell's equations is: 



Ei,± = 



E, 



i,±e 



iki,± -r—itot 



(Bl) 
(B2) 



where + represents a wave traveling in the positive z 
direction at some otherwise arbitrary angle, — represents 
a wave traveling in the negative z direction, and i is the 
layer number. Then in CGS: 



Hi,± = J— k h± x Ei- 



— ki.± x Ei.± 



CGS (B3) 



where ki y ± is a unit vector parallel to the propagation 
direction, ej is the complex dielectric constant of the zth 
layer and includes all conductivity effects and Zi = 

is the relative impedance of the medium and like e may 
be complex. 

Now that we have a solution inside a layer we need 
to match solutions for different layers at the boundaries 
between layers. The curl E Maxwell equation provides 
one boundary condition: 



Ei 



x z 



E, 



+i x z. 



(B4) 
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The curl H Maxwell equation provides the other bound- 
ary condition: 



Hi x z = H, 



+1 x z. 



(B5) 



To apply the foregoing consider first S polarization (E 
perpendicular to the plane of incidence) at each side of a 
boundary In such an instance the E field is tangent to 
the boundary and if we orient the x and y axes so that 
E is also perpendicular to the y axis, then, from Eq. B4: 



£i, + ,x ~\~ £i. 



£i+l,+,x + £i+l, — ,x- 



(B6) 



Likewise, substituting Eq. |B3| into Eq. B5 produces: 
cos(0j) 



-(£ 



Si. 



w+i 



-{£ 



i+l, + .x 



£, 



i+l,—,x ) 



(B7) 

where 6i is the angle of incidence measured of the z axis. 
Defining the incident impedance for S polarization as: 



cos(6»i) 



(B8) 



Eq. B7 becomes 

Zs,i(£i, + ,x — £i,-,x) = Zs : i+l(£i+l, + : x — £i+l,-,x)- (B9) 

and 



Combining Eqs. B6 



£i+l,+ 
£%+!,- 



|B9| into a matrix form: 



£,- 



where 



2Z, 



Z s ,i+1 Z s 
Z s ,i+1 Z s 



J s,+i 
J s,+i 



(BIO) 



(Bll) 



For P polarization (E field || to plane of incidence) 
define the incident impedance as: 



then following the same procedure: 



P, 



(B13) 



where 



P 



1 I ^p,i+l + Zp, : 



IZr, 



z 



z. 



p,i+l 



Z>n 



p,i+i 



Z n 



(B14) 



The E field can be determined from H using Eq. BS 



Equations |BJ, |B1C| , and §11| and Eqs. §1J, |B13 
and B14 provide the incident and reflected amplitudes 
in the i + 1 layer at the interface in terms of the ampli- 
tudes in the i layer at the interface for S and P polar- 
ization respectively. It remains to develop the equations 
to propagate the amplitudes across lay er i of thickness 
di. This is already provided by Eqs. Bl and B2 which in 
matrix form become: 



£i.+ {z = di) 
£i-{z = di) 







^0 



£ t .+ (z = 0) 
£,-,-(* = 0) 



= Ui 



£i,+( z = °) 

£ it -(z = 0) 



So, as an example, given a set of layers numbered 1 
to n from left to right the equation relating the incident 
and reflected amplitudes on left side of the set to the 
right side for S polarization is: 



Sii-i.iiU n _iS n _2,ti-i ■ ■ • Ui + iS 



£i,+( le fty 

" 1 1 £i. (left), 

£ rh+ (right)\ 
£ n - (right) J 



cos(9 l )Z i 



(B12) 
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FIG. 1: Overall schematic of the optical path. The the dotted line represents the C02laser beam. The CO2 laser and the first 
four optic elements are actually located on a separate mount detached from the optical table. 




FIG. 2: Wedged sample and beam path. The paths of the 
reflected beams are sufficiently separated by the 2° wedge to 
prevent their entering the detection system or reentering the 
laser. 



4 



14 



Direction of 
Laser Beam 



Steady Pin 



Indexing Hole 




Heater 



Nylon Screws 



BaF Wedged Window 
Background Sample 



Flexible 
Thermal Bond 




Copper Sample 
Carrier 



BSCCO Insulating Spacers 

Sample Film Retaining Spring 



FIG. 3: Sample mount assembly 
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FIG. 4: The 2tj m and 3uj m normalized signals for a Cu film as a function of magnetic field at 949 cm 1 and 210 K. 
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FIG. 5: The Re[#/] as a function of frequency squared at 8 T for a GaAs sample. The solid line represents a uj 2 -&t. The 
semi-insulating sample is wedged with an average thickness of 0.42 mm. 



